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This  report  summarizes  the  results  obtained  during  an 
investigation  directed  toward  the  use  of  adpative  antenna 
array  beamforming  techniques  in  HF  backscatter  radar  systems. 

The  work  reported  covers  the  period  January  1,  1973  through 
February  28,  1974  and  is  a  continuation  of  previously  reported 
studies.  Data  for  this  investigation  were  obtained  using  the  Side- 
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Los  Banos  in  the  central  valley  of  California  and  is 
operated  by  Stanford  Research  Institute,  Menlo  Park,  Cali¬ 
fornia  Simultaneous  digitized  recordings  were  taken  at 
our  stibrray  outputs,  with  each  subarray  consisting  of  32 
vertical  monopole  elements  which  were  combined  using  a 
conventional  time  delay  and  sum  beamforming  system.  Total 
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Ihree  transmit /receive  formats  were  used  during  these 
investigations:  l)  CIV  signals  transmitted  at  Bearden, 

;™nSKS  at  3  frecluenc>’  selected  to  provide  one-hop 
onosphenc  propagation  were  received  at  Los  Banos  in  the 
1  resence  of  strong  teletype  interference.  2)  CW  signals 
transmitted  at  Lost  Hills,  California  were  used  to  provide 
one-hop  backscatter  geometry.  A  repeater,  located  at 
i earden,  provided  100  Hz  DSB  modulation  for  purposes  of 
arget  simulation.  3)  FM/CIV  radar  format  signals  were 
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ity  of  Colorado  using  a  CDC  6400  computer.  No  real-time 
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Comparison  of  adaptive  and  conventional  beamforming  tech¬ 
niques  on  these  data  demonstrated  that  an  output  signal-to- 
noxse  ratio  (SNR)  .improvement  of  25-  30  dB  mayPbe  achieved 
using  adaptive  techniques  m  the  presence  of  CIV  transmissions 
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operate  in  real-time  on  the  SRI  IVARF  array  and  that  such 
metnods  will  significantly  improve  overall  radar  performance 


16b 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGEflfh.n  Dmtm  En,ertd) 


ADAPTIVE  ARRAY  PROCESSING  Op  HF 
BACKSCATTER  RADAR  SIGNALS 

Lloyd  J.  Griffiths 


Contractor:  University  of  Colorado 
Contract  Number:  F306Q2-72-C-0386 
Effective  Date  of  Contract:  1  January  1972 
Contract  Expiration  Date:  31  July  19/4 

Amount  of  Contract:  $52/i30.00 
Program  Code  Number:  4E20 


Principal  Investigator: 

Phone: 


Prof  L.  J.  Griffiths 
303  443-2211,  ext  7653 


Project  Engineer: 

Phone: 


Leonard  Strauss 
315  330-3055 


Approved  for  public  release; 
distribution  unlimited. 


This  research  was  supported  by  the 
Defense  Advanced  Research  Projects 
Agency  of  the  Department  of  Defense 
and  was  monitored  by  Leonard  Strauss 
RADC  (0CSE)  GAFB  NY  13441  under 
Contract  F30602-72-C-0386. 


>b 


ABSTRACT 


This  report  summarizes  the  results  obtained  during  an 
investigation  directed  toward  the  use  of  adaptive  antenna 
array  beamforning  techniques  in  HF  backscatter  radar  systems. 
The  work  reported  covers  the  period  January  1,  1973  through 
February  2S,  1974  and  is  a  continuation  of  previously 
reported  studies. 

Data  for  this  investigation  were  obtained  using  the 
Wide-Aperture  !1F  Radio  Research  Facility  which  is  located 
near  Los  Banos  in  the  central  valley  of  California  and  is 
operated  by  Stanford  Research  Institute,  Menlo  Park,  Cali¬ 
fornia.  Simultaneous  digitized  recordings  were  taken  at 
four  subarray  outputs,  with  each  subarray  consisting  of 
32  vertical  monopole  elements  which  were  combined  using  a 
conventional  time  delay  and  sum  beamforming  system.  Total 
aperture  for  the  study  was  1.3  km  Individual  receivers 
with  1  kHz  bandwidth  and  common  local  oscillator  were 
used  at  the  four  subarray  outputs.  Transmitting  facilities 
at  Lost  Hills,  California  and  Bearden,  Arkansas  were  used 
during  the  experiment. 

Three  transmit/reccive  formats  were  used  during  these 
investigations:  1)  CW  signals  transmitted  at  Bearden,  Arkan¬ 
sas  at  a  frequency  selected  to  provide  one~hop  ionospheric 
propagation  wer^  received  at  Los  Banos  in  the  presence  of 
strong  teletype  interference.  2)  CW  signals  transmitted 
transmitted  at  Lost  Hills,  California  were  used  to  provide 
one-hop  backscatter  geometry.  A  repeater,  located  at  Bearden, 
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provided  100  Hz  DSB  modulation  for  purposes  of  target  sim¬ 
ulation.  3)  FM/CW  radar  format  signals  were  transmitted 
from  Lost  Hills  and  the  radar  was  used  to  search  for  targets 
of  opportunity. 

In  all  three  cases,  data  were  processed  at  the  Univer¬ 
sity  of  Colorado  using  a  CDC  6400  computer.  No  real-time 
processing  was  attempted  during  the  course  of  this  study. 
Comparison  of  adaptive  and  conventional  beamforming  tech¬ 
niques  on  these  data  demonstrated  that  an  output  signal-to- 
noise  ratio  (SNR)  improvement  of  25-30  dB  may  be  achieved 
using  adaptive  techniques  in  the  presence  of  0:1  transmissions 
This  improvement  was  achieved  for  both  the  forward  scatter 
and  backseatter  Geometries.  Adaptation  performance  in  the 
presence  of  clutter  was  identical  to  that  in  an  interfer¬ 
ence  only  environment.  Total  time  for  adaptation  from  zero 
coefficient  values  to  steady-state  optimum  values  was  shown 
to  be  less  than  40  msec.  Substantial  noise  reduction  im¬ 
provement  in  the  range/doppler  display  obtained  from  an 
1  M/0<1  transmitted  signal  was  demonstrated. 

The  conclusion  drawn  from  this  study  is  that  simple 
adaptive  beanforming  methods  can  easily  be  employed  to 
operate  in  real-time  on  the  SRI  WARF  array  and  that  such 
metnods  will  significantly  improve  overall  radar  performance. 
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I.  SUMMARY 


This  report  outlines  the  technical  results  obtained 
under  RADC  Contract  No.  F  3060 2-72-C-O 38C  during  the  period 
January  1,  1973  to  February  28,  1974.  The  work  reported  is 
a  continuation  of  that  previous  reported.  ^  This  section 
presents  a  summary  of  the  primary  results  obtained  during 
the  course  of  this  study  and  the  conclusions  drawn  from 
these  results.  Detailed  information  is  presented  in  the 
sections  following. 

The  objective  of  this  work  was  to  obtain  quantitative 
measurements  of  the  improvement  provided  by  an  adaptive  array 
beamforming  system  in  HF  radar  applications.  Results  were 
obtained  using  experimental  data  recorded  in  a  relatively 
noisy  HF  environment.  This  work  is  a  continuation  of  a  pre¬ 
viously  reported  study Il]  which  utilized  CW,  one  hop,  for¬ 
ward  scatter  ionospheric  signals,  and  extends  these  results 
to  the  case  of  realistic  HF  backscatter  radar  environment. 

Data  j.or  this  investigation  were  recorded  using  the  Wiue- 
Aperture  HF  Radio  Research  Facility  (WARF)  operated  by 
Stanford  Research  Institute,  Menlo  Park,  California.  A 
complete  description  of  this  facility  is  given  in  Reference  2. 
An  improved  receiving  system  was  employed  for  the  present 
investigation.  This  system  had  1  kHz  bandwidth  receivers  at 
each  of  four  subarray  outputs  at  the  WARF  receiving  site  at 

1 


Los  Banos . 


Simultaneous  sample  and  hold  circuitry  was  in¬ 
stalled  to  provide  a  2  kHz  sampling  rate  at  each  receiver 
output.  All  data  were  recorded  on  seven  track  computer 
compatible  magnetic  tape  and  nailed  to  the  University  of 
Colorado  for  processing.  Computations  were  carried  out 
using  the  CDC  6400  facility  at  Colorado.  Data  for  the  pre¬ 
sent  study  were  recorded  during  August,  1973. 

Data  were  taken  simultaneously  from  four  adjacent  sub- 
irray  outputs  of  the  WARF  receiving  array.  Each  subarray 
contains  32  vertical  monopole  elements  at  10  meter  spacing. 
These  elements  are  combined  using  a  remotely  controlled  time 
delay  and  sum  beamforming  system.  The  total  aperture 
spanned  by  the  four  subarrays  is  1.3  km.  Transmitting 
facilities  at  Lost  Hills,  California  were  used  to  generate 
the  radar  signals  and  a  low  power  repeater  at  3earden,  Arkan¬ 
sas  was  used  for  purposes  of  target  simulation.  in  all 
tests,  transmitter  frequencies  were  selected  so  as  to  pro¬ 
vide  a  one-hop  ionospheric  path  between  Bearden  and  Los 
Banos.  All  tests  were  conducted  with  the  subarray  elements 
steered  in  the  boresight  --  i.e.,  90  degrees  true  -- 
direction.  Transmitter  frequencies  were  intentionally  cho¬ 
sen  near  high  interference  regions  of  the  HF  spectrum  for 
purposes  of  determining  the  SNR  improvement  capabilities  of 
the  adaptive  beamforming  system.  All  tests  were  conducted 
at  approximately  local  noon  at  the  ionospheric  midpath  point. 

Three  different  experimental  configurations  were  imple¬ 
mented  during  the  course  of  this  study: 
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1*  forward  scatter  test— low  power  CW  signals  were 
transmitted  from  Bearden  and  received  at  Los  Banos. 
Frequencies  were  selected  to  provide  both  high  and  low 
SNR  conditions. 

2.  CW  backscatter  test- fixed  frequency  transmissions 
were  generated  at  Lost  Hills  and  received  in  a  one- 
hop  backscatter  mode  at  Los  Banos.  A  repeater  at 
Bearden  was  used  to  produce  100  Hz  sideband  modulation 
for  purposes  of  simulating  a  target.  Both  high  and  low 
SNR  conditions  were  provided. 

3.  FM/CW  backscatter  test  -  a  radar  waveform  format 
consisting  of  a  sweep  over  8  1/3  kHz  at  a  60  Hz  repi- 
tition  rate  was  transmitted  from  Lost  Hills.  The  Los 
Banos  array  was  then  used  to  locate  targets  of  oppor¬ 
tunity  in  the  backscatter  signal.  Once  such  a  target 
had  been  located,  data  were  recorded  at  the  eight  sub¬ 
array  outputs. 

These  data  were  processed  at  Colorado  to  evaluate  the  per¬ 
formance  improvement  provided  by  an  adaptive  beamforming  sys¬ 
tem.  In  each  case,  an  adapted  beam  output  was  computed  at 
a  given  data  sample  as  the  sum  of  delayed  and  weighted  sig¬ 
nals  out  of  the  four  subarray  receivers.  The  number  of  de¬ 
layed  signals  from  each  subarray  —  i.e.,  the  total  memory 
duration  per  subarray  —  was  preselected  as  1,3,5,  or  9 
samples.  Thus,  the  total  number  of  variable  weights  in  the 
adaptive  array  processor  was  4,12,20,  or  36.  These  weights 
were  updated  as  each  new  data  sample  was  received  using 
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the  P-vector  adaptive  algorithm  described  in  the  sections 
following.  The  resulting  beam  output  was  then  compared  with 
a  conventional  beam  output  formed  by  applyinq  fixed  Tchebsychev 
weights  to  the  four  subarny  outputs  and  summing  the  result. 

Average  output  SNR  improvement  of  adaptive  over  conven¬ 
tional  beamforming  was  measured  for  both  the  forward  and  back- 
scatter  CW  tests  described  above.  The  results  obtained  were 
as  follows: 


1.  CW  forward  scatter  at 
subarray  output) 

no.  of  wts/subarray 
1 
3 
5 
9 

2.  CW  backscatter  at  high 

5 

3.  CW  backscatter  at  low 

5 


received  SNR  (-25  dB  at  a 

SNR  improvement 
3.2  dB 

5. 5  dB 

9.6  dB 
11.1  dB 

received  SNR  (-38  dB) 

25.5  d3 
eceived  SNR  (-44  dB) 

27.9  dB 


It  should  be  noted  that  the  received  SNR's  shown  above  for 
the  backscatter  signals  are  computed  as  signal-to-noisc  plus 
interference  plus  clutter  ratios.  The  total  time  to  reach 
steady-state  adaptation  was  measured  for  each  of  the  above 
cases.  This  measurement  was  carried  out  by  setting  the 
adaptive  coefficients  to  zero  and  monitoring  the  output  SNR 
as  a  function  of  time  as  the  weights  adapted.  The  results 
varied  between  15  m  sec  for  9  weights  per  subarray  in  the 


4 


fo ’-ward-scatter  experiment  to  30  m  sec  xor  5  weights  in  the 
low  SNR  backscatter  test. 


Comparisons  were  also  derived  for  the  FM/CW  tests. 

For  these  data,  range/doppler  maps  were  computed  using  a 
0.5  sec.  integration  interval.  Maps  were  computed  for  the 
sane  data  set  using  both  fixed,  conventional  beamforming  at 
the  four  subarray  outputs  and  a  time  variable,  continuously 
adapting  beamforner.  In  all  cases,  the  adaptive  range/ 
doppler  displays  exhibit  substantial  noise  improvement. 
Specific  ox amnios  are  » resented  in  the  sections  following. 

Several  conclusions  may  be  drawn  from  the  results  of 
this  study: 

1.  Output  SNR  improvements  of  20-30  d.3  may  be  achieved 
using  four  adaptive  subarrays. 

2.  Total  time  to  adapt  is  less  than  40  m  sec. 

3.  Five  adaptive  weights  per  subarray  is  sufficient  to 
provide  substantial  improvement. 

4.  Adaptive  improvement  against  clutter  is  identical 
to  that  achieved  for  interference  --  i.e.,  the  fact 
that  the  clutter  is  correlated  with  the  desired  signal 
does  not  degrade  adaptive  performance. 

5.  The  P-vector  algorithm  can  be  used  effectively  for 
FM/CW  radar  formats. 

Thus,  adaptive  processing  of  HF  radar  array  signals  appears 
to  offer  considerable  improvement  over  conventional  beam¬ 
forming  systems,  Further  investigations  should  now  be  car¬ 
ried  out  to  implement  and  demonstrate  an  on-line,  real-time 
adaptive  system. 
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II.  INTRODUCTION 


The  investigation  described  herein  is  a  continuation  of 
a  program  which  was  initiated  early  in  1972.  Previous  re¬ 
sults  obtained  during  this  program  are  summarized  in  Refer¬ 
ences  1  and  3.  These  reports  contain  bibliographies  and 
background  information  relevant  to  the  general  area  of 
adapt- ive  beamforming  techniques.  The  primary  objective  of 
the  overall  program  is  to  define  the  improvement  provided  by 
an  adaptive  beamforming  system  in  HF  radar  applications.  The 
work  has  been  directed  toward  an  investigation  of  the  use 
of  adaptive  methods  in  the  radar  receiving  array  only.  Thus, 
the  transmitting  system  is  presumed  to  consist  of  both  fixed 
antenna  hardware  and  fixed  transmitter  waveform  format.  In 
effect,  the  proposed  adaptive  array  beamformer  is  to  be 
considered  as  a  component  which  will  be  added  to  an  existing 
radar  system  with  minimal  change  in  the  overall  configuration . 
The  questions  which  must  be  answered  during  the  course  of 
this  program  are  1)  What  is  the  improvement  offered  by  the 
proposed  system  and  2)  What  are  the  costs  associated  with 
its  implementation. 

The  specific  objective  of  the  research  reported  in 
this  report  was  to  obtain  quantitative  measurements  which 
provide  an  answer  to  the  first  of  these  questions.  The 
measurements  were  taken  under  conditions  representative  of 


tho  environment  in  which  an  HF  system  must  operate.  Spe¬ 
cifically,  one-hop  backscatter  signals  generated  by  an 
FM/CW  transmitter  format  were  received  in  the  presence  of 
relatively  hiah  noise  and  interference  levels.  Adaptive 
beam forming  techniques  are  ideally  suited  to  these  conditions 
and  can  provide  time  variable  null-tracking  of  both  the  fre¬ 
quency  and  spatial  direction  of  arrival  of  interference  sour¬ 
ces. 

A  partial  answer  to  the  second  question  stated  above  may 
also  be  inferred  from  the  results  of  the  present  study.  The 
specific  equipment  configuration  used  to  obtain  the  results 
consisted  of  an  all  digital  system  similar  to  that  shown  in 
Fin.  1.  '[’he  cost  of  this  system  is  easily  ascertained  and 

bears  a  direct  relationship  to  the  number  of  adaDtive  sub¬ 
arrays  used  and  the  complexity  of  the  adaptive  processing 
algorithm.  Further  investigations,  however,  are  required  to 
evaluate  hybrid  configurations  of  the  type  shown  in  Fia.  2 
which  may  have  significant  advantages  for  wide  bandwidth 
and/or  multiple  element  systems. 


recv.  k-  ■  4  ADC 
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Fig.  1  General  form  of  all  digital  adaptive  receivina 
array  beamforming  system. 


III.  EXPERIMENTAL  FACILITIES  AND  DATA  BASE 


A.  EQUIPMENT 

All  data  for  the  present  study  were  recorded  at  the 
Los  Banos  receiving  site  of  the  WARF  facility  operated  by 
the  Ionospheric  Dynamics  Group  at  Stanford  Research  Insti¬ 
tute,  Menlo  Park,  California.  Active  transmitters  at 
Bearden,  Arkansas  and  Lost  Hills,  California  were  used 
during  the  course  of  the  study.  Figure  3  illustrates  the 
relative  locations  of  these  facilities  and  summarizes  their 
parameters.  All  tests  were  conducted  with  both  the  Los 
Banos  and  the  Lost  Hills  arrays  steered  in  the  90°  true 
direction. 

1  tie  receiving  system  at  Los  Banos  utilized  the  outputs 
of  four  spatially  separated  adjacent  subarrays.  Each  output 
consisted  of  the  sum  of  32  vertical  monopole  elements  with  a 
10  meter  spacina  between  elements.  The  array  is  linear  and 
oriented  in  the  North-South  direction.  Thus,  a  given  sub¬ 
array  aperture  was  320  meters  and  the  total  aperture  spanned 
was  1.28  km.  It  should  be  noted  that  this  aperture  represents 
one  half  of  the  total  available  at  Los  Banos  and  that  the 
entire  system  contains  eight  subarray  outputs.  Equipmental 
difficulties  at  the  time  of  the  tests,  however,  prevented 
utilization  of  the  full  array. 

figure  4  illustrates  the  receiving  and  recording  system 
used  m  this  work.  The  subarray  output  signals  were  de- 
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FIGURE  3  Location  map. 


W1DE-APERTURF  RECEIVING  STATION  AT  LOS  BANOS.  CALIFORNIA: 

256  vortical  monopolo  antenna  elements  in  8  sub-arrays. 

Total  length  of  array,  2  5  ktn. 

Frequency  range.  0  to  27  MHz 
Bcamwidth  at  20  MHz.  0  5  deg. 

Array  is  bidirectional,  broadside 
Boresight.  90  deg  and  270  deg  true 

Direction  of  main  lobe  of  sensitivity  pattern  may  be  repetitively  slewed  ±  16  deg  in 
0.25  deg  steps 
Terrain  eastward.  Oat. 

TRANSMITTING  AND  REPEATING  STATION  AT  BEARDEN.  ARKANSAS 
Average  power  capability.  I  kW 
Antennas 

One  single-element  ITT  FTM: 

Two  horizontally  polarized  Collins  LPA's  (rotatable:  13  dB  gain). 

TRANSMITTING  STATION  AT  LOST  HILLS.  CALIFORNIA: 

Average  power  capability.  30  kW. 

Antenna.  18-element  ITT  LDAA. 

Normal  gam  18  dB 
Midband  bcvinwidth.  6  deg. 

Boresight,  00  deg  true,  viewable  t  32  deg  in  4  deg  steps 


11 


modulated  using  four  identical  mixers  and  receivers  driven 
by  a  common  local  oscillator.  The  receivers  used  were  in¬ 
stalled  in  March,  1973  and  contained  several  improvements 
over  those  used  for  the  previous  study  A  receiver  band¬ 

width  of  1.0  kH*  was  selected.  The  detected  baseband  output 
signals  served  as  inputs  to  simultaneous  sample  and  hold  cir¬ 
cuits  which  were  then  digitized  using  a  multiplexed  12  bit 
analog  to  digital  converter  and  stored  on  magnetic  tape.  Thus, 
simultaneous  samples  were  recorded  from  each  of  the  four  sub¬ 
arrays  and  the  skew  correction  used  previously^  to  compensate 
for  a  successive  element  A/D  conversion  was  not  required.  The 
effective  sampling  rate  of  the  system  was  2.0  kHz  per  subarray 
output.  The  magnetic  tapes  produced  in  this  manner  were  later 
converted  to  CDC  6400  compatible  format  and  shipped  to  the 
University  of  Colorado  for  further  processing. 

B .  DATA  BASE 

The  magnetic  tape  format  used  to  collect  data  was  seg¬ 
mented  into  individual  files  with  each  file  representing  a 
time  continuous  recording  c.  f  a  particular  experimental  con¬ 
figuration.  A  given  data  file  contained  between  2  and  1000 
physical  magnetic  tape  records  with  each  record  consisting 
of  64  milliseconds  of  data  from  the  four  subarray  outputs. 

Table  I  is  the  data  log  written  on  August  17,  1974. 

Local  time  for  the  tests  was  12:00  to  15:00  PST.  Twelve 
files  of  data  were  recorded  on  this  date  and  may  be  grouped 
into  the  following  four  experimental  configurations: 


13 


during  12:00  -  15:00  PST 


test 


1.  File  1  Calibration: 

Simultaneous,  in  phase,  equal  amplitude  sinusoidal  test  sig¬ 
nals  were  inserted  into  the  system  shown  in  Fig.  4  at  the 
locations  normally  used  to  connect  the  subarray  output  signals 
to  the  double  balanced  mixers.  These  signals  are  used  for  in¬ 
dividual  gain/phase  matching  of  the  receivers  and  serve  as  an 
indication  of  the  system  spurious  levels. 

2.  Files  2-6  Forward  scatter  CW  tests: 

Transmitting  facilities  located  at  Bearden  were  used  to  gen¬ 
erate  signals  at  17.124  MHz.  This  frequency  was  selected  to 
provide  one-hop  ionospheric  propagation  over  the  Bearden  to  Los 
Banos  path  as  well  as  a  moderate  amount  of  interference  in  the 
1  kHz  operating  bandwidth.  The  data  recorded  on  files  3  and  4 
consisted  of  CW  transmissions  with  the  carrier  offset  by  500 
HZ  so  as  to  appear  in  the  center  of  the  1000  Hz  analysis  band. 
Signal- to-noise  ratio  was  varied  by  reducing  the  transmitter 
power.  Files  5  and  6  were  recorded  in  a  similar  fashion  with 
the  addition  of  a  1000  Hz  broadband  noise  modulation  imposed 
on  the  CW  carrier  for  purposes  of  simulating  a  spread  spectrum 
source.  The  objective  of  this  series  of  experiments  was  to 
obtain  data  under  conditions  similar  to  those  taken  during 
the  1972  study  and  to  use  the  data  as  a  basis  for  comparison. 

3.  Files  7-8  Backscatter  CW  tests: 

The  Lost  Hills  facility  was  used  to  transmit  CW  signals  at 
16.295  MHz-  A  low  power  repeater  at  Bearden  induced  a  +100  Hz 
offet  prior  to  retransmission  of  the  received  signal.  The 
operating  frequency  provided  one-hop  propagation  over  the 
path.  Received  SNR  was  varied  during  this  experiment  by  re- 
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ducing  the  repeater  output  power.  These  data  were  collected 
for  the  purposes  of  comparing  adaptive  eprformance  in  the 
presence  of  clutter  with  that  obtained  under  conditions  of 
interference  only  (Files  2-6). 

4.  Files  9-12  FM/CW  radar  experiments: 

The  purpose  of  these  experiments  was  to  collect  data  using 
transmitter  waveforms  identical  to  those  which  would  be  used 
in  an  operational  HF  radar  system.  The  waveform  transmitted 
from  Lost  Hills  consisted  of  a  sawtooth  linear  FM  sweep  over 
a  8  1/3  kHz  frequency  range,  centered  at  16.3  MHz,  and  a 
sweep  repetition  rate  of  60  Hz  as  shown  in  Fig  5.  A  repeater 
at  Bearden  was  employed  to  simulate  a  15  Hz.  doppler  offset 
target.  The  experiments  consisted  of  tests  under  both  low  and 
high  SNR  conditions  and  a  "pseudo-operational"  experiment  at 
11.64  MHz  with  the  repeater  turned  off.  During  this  test, 
the  radar  was  used  to  search  for  targets  of  opportunity.  As 
shown  in  Table  I,  the  sampling  rate  was  reduced  to  1.920  kHz 
to  provide  synchronization  with  the  60  Hz  chirp  repetition 
rate . 

C.  DATA  PROCESSING  AND  EVALUATION 

One  method  which  may  be  used  to  examine  spatially  dis¬ 
tributed  digital  data  which  has  been  obtained  from  uniformly 
spaced  linear  array  is  the  two-dimensional  digital  Fourier 
transform  (2D  DFT)  described  in  a  previous  report^.  In  this 
technique,  a  two-dimensional  data  matrix  A  is  formed  as 


follows : 


^(k)  x^(k-l)  ...  x^(k«-N+l) 

2(k)  x2(k-l)  ...  x2(k-N+l) 

3  (k) 

4(k)  x^  (k-N+1) 

t  h 

where  x^  (k-j)  =  sample  at  time  (k-j)  from  x  subarray  and  N  is 

the  number  of  samples  per  element  processed.  Transforming  the 

rows  and  columns  of  A  produces  a  new  matrix  in  which  the 

t  hi 

elements  of  the  %  ‘  row  may  be  interpreted  as  the  frequency 

t  hi 

spectrum  for  signals  received  from  the  P.  resolvable  azimuth. 
In  a  four  element  array,  there  are  four  resolvable  azimuths 
spaced  by 


where 

C  =  velocity  of  propagation 
L  =  no.  of  elements  (subarrays) 
f  -  received  frequency 
d  =  element  spacing. 

For  the  4  subarray  system  used  in  the  present  study  the 

spacing  at  17  MHz  is  0.80  degrees. 

Figure  6  shows  the  results  obtained  with  this  technique 

for  the  calibration  signal  recorded  in  File  1.  A  total  of  256 

samples  per  element  were  processed  to  produce  the  display, 

thus  providing  a  frequency  resolution  of  about  8  Hz.  The 

first  eight  zimuthal  lines  shown  were  obtained  by  interpolation 

th 

of  the  four  lines  provided  by  the  data.  (The  91"  line  is  a 
reproduction  of  the  first  included  to  maintain  symmetry.)  A 
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FREQUENCY 


Fig. 


a)  Display  gain  =1.0 


b)  Display  gain  =  100.0 


Azimuth/ frequency  display  for  calibrate  signal. 
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nine  line  display  was  generated  for  purposes  of  comparison 
with,  the  previous  work  which  resulted  from  eight  element 

cl  cl  1 3  • 

Fig.  6a)  has  been  normalized  to  the  maximum  computed 
value,  and  6b)  is  the  same  display  with  all  values  multiplied 
by  a  factor  of  100.  Large  amplitude  values  in  Fig.  6b)  have 
been  clipped.  The  calibration  signal  is  located  at  500  Hz  and 
has  an  apparent  broadside  direction  of  arrival.  The  sidelobes 
shown  are  produced  by  a  raised  cosine  weighting  which  was 
applied  to  the  data  prior  to  2D-DFT .  Note  that  relatively  few 
spurious  signals  are  present  in  the  data.  The  noise  floor 
exhibits  the  same  spatial  characteristics  as  the  calibration 
signal  and  was  most  likely  produced  by  the  calibration  signal 
generator.  These  results  represent  a  considerable  improve¬ 
ment  over  "he  receiving  system  used  in  the  previous  study.  The 
signal  amplitude  relative  phase,  ana  D.C.  offset  measured  at 
each  receiver  output  during  calibration  is  summarized  in 
Table  II. 

Table  II :  Calibration  Summary 

Phase  DC  Offset 

0  0.059 

0.021°  0.015 

-0.042°  0.012 

0.005°  0.079 

Figs.  7-10  illustrate  2D  DFT ' s  of  some  typical  sample  of 
the  raw  data  recorded  during  the  present  study.  The  forward- 


channel 

1 

2 

3 

4 


Amplitude 
4.88 
4.86 
4 . 84 
4  .  87 
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scatter  CW  transforms  shown  in  Fig.  7  were  computed  using  a 
total  of  512  time  samples  from  each  subarray  output.  The 
resulting  frequency  resolution  in  the  display  is  4  Hz  and 
enables  identification  of  the  desired  signal  in  the  high  SNR 
case.  Similar  labelling  for  the  low  SNR  case,  however,  is  not 
readily  obtained.  Many  of  the  remaining  transforms  presented 
in  this  report  were  computed  using  only  128  samples  per  ele¬ 
ment  and  a  resulting  resolution  of  16  Hz.  The  transform 
equivalent  to  Fig.  7b)  is  presented  in  Fig.  8  for  comparison. 
Note  that  although  much  of  the  fine  scale  frequency  structure 
can  no  longer  be  identified,  spatial  and  frequency  domain 
characteristics  of  the  interference  field  are  readily  apparent. 
The  broadband  signal  shown  in  the  600-750  Hz  frequency  range 
was  a  teletype  interference  identified  at  the  time  date  were 
recorded. 

Figs.  9  and  10  illustrate  the  clutter-dominated  data  ob¬ 
served  during  the  CW  backscatter  tests.  The  repeater  offset 
signal  can  be  seen  in  the  high  SNR  experiment  (arrow)  and 
interference  lines  appear  in  the  low  SNR  display.  Transform 
results  for  the  FM/CW  experiment  are  shown  in  Fig.  11.  The 
multiplicity  of  clutter  lines  which  appear  at  60  Hz  separa¬ 
tions  corresponds  to  the  chirp  repetition  rate. 
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Fig.  7  a) 

Azimuth /frequency 
display  for  high 
SNR  forward  scatter 
data. 


Fig.  7  b) 

Azimuth/ frequency 
display  for  low 
SNR  forward  scatter 
data. 


Fig.  8 

Azimuth/ frequency 
display  for  low 
SNR  forward  scatter 
data  using  128  pt. 
trans  form. 
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a)  Di.  -'lay  gain  =  1.0 


b)  Display  gain  =  3.0 


Azimuth/ frequency  display  for  high  SNR  backscatter  data 


FREQUENCY 


a)  Display  gain  =  1.0 


Fig.  10  Azimuth/frequoncy  display  for  low  SNR  backscatter  data. 


24 


IV.  ARRAY  PROCESSING  THEORY 

The  basic  theory  of  array  processing  of  interest  in  the 
present  study  has  been  developed  in  a  previous  report . 

The  theory  presented  below  extends  this  work  to  include  a  cal¬ 
culation  of  the  adaptive  convergence  time  and  the  effects  of 
adaptation  in  the  presence  of  an  FM/CW  transmitted  waveform. 

A  GAIN-PHASE  BEAMFORMING  VS.  FREQUENCY-DEPENDENT  BEAMFORMING 

The  beamf orming  system  presented  in  this  study  is  an 
all  digital  system  in  which  successive  samples  from  the  down- 
converted  subarray  outputs  are  weighted  and  summed  to  pro¬ 
duce  a  single  sampled  beam  output  signal.  Figure  12  illustrates 
the  basic  processor.  Note  that  each  antenna  "element"  is  ac¬ 
tually  a  32  element  subarray.  This  processor  is  capable  of 
filtering  signals  incident  on  the  array  both  in  terms  of  the 
frequency  of  the  received  signal  and  the  direction  of  arrival. 
Thus,  the  beamformer  acts  as  a  combined  spatial/temporal  filter. 

To  illustrate,  assume  a  planar  sinusoidal  signal  at  fre¬ 
quency  fc+fm  arrives  from  azimuthal  angle  0.  Let  the  rela¬ 
tive  phase  delay,  in  radians,  between  two  successive  elements 
be  denoted  by  A0.  We  further  assume  that  the  mixers  are 
driven  with  a  common  local  oscillator  at  frequency  Hz  and 
that  fc  is  within  the  operating  bandwidth  of  the  system 
(i.e.,  '  f  jp  |  £  500  Hz).  The  signal  observed  after  down  con- 
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12  Tapped- del ay- line  digital  array  processing  configuration. 


version  is  then  a  sinusoidal  signal  at  frequency  f  and  the 
relative  phase  between  successive  outputs  remains  at  AG 
radians  due  to  the  use  of  a  common  local  oscillator.  Using 
phasar  notation  the  k^1  sample  out  of  the  ifc^  element  x^  (k) 
may  be  represented  as 


Equation  (6)  illustrates  the  spatial-temporal  duality 


of  the  digital  beamformer  shown  in  Fig.  12.  Complete 
symmetry  is  obtained  for  the  case  L  =  K  — -  i.e.  ,  the  number 
of  digital  coefficients  per  element  is  equal  to  the  number 
of  spatial  elements.  To  illustrate  the  duality,  assume  L  =  K 
and  that  a  set  of  coefficients  (wi  £}  has  been  determined 
such  that  the  array  frequency  response  to  signals  incident 
from  direction  0^  and  having  a  H/2  phase  shift  between 
successive  elements  is  as  shown  below. 


Fig.  13  Assumed  processor  frequency  response  for  signal 

arriving  from  direction  0  . 

o 

Further  assumed  that  a  new  processor  is  implemented  with  com¬ 
plementary  coefficients  {5^=  w^.].  This  system  will  then 
have  a  spatial  beam  pattern  for  signals  incident  at  frequency 
%AT  =  n/2  which  is  identical  to  that  shown  in  Fig.  13. 
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Fig.  14  Processor  spatial  beam  pattern  for  system  con¬ 
taining  complementary  coefficients. 


Note  that  the  horizontal  axis  in  this  plot  is  the  inter¬ 
element  progressive  phase  shift  and  is  not  necessarily  the 
direction  of  arrival. 

An  alternative  all-digital  array  processor  which  uti¬ 
lizes  quadrature  receivers  is  shown  in  Fig.  15.  For  the 
narrow  bandwidths  of  interest  in  the  present  study,  it  is 


presumed  that  the  quadrature  relationship  is  preserved  for  all 

i 

received  signals.  Adjustment  of  the  coefficients  in  this 
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system  will  provide  only  spatial  processing  —  i.e.  ,  frequency 

selectivity  within  the  bandwidth  cannot  be  achieved. 

To  illustrate,  using  the  assumptions  given  prior  to  (3) 
t  h 

above,  the  k  sample  of  the  downcon verted  in-phase  and  quad¬ 
rature  components  p^k)  and  q^  ( k)  ,  respectively,  observed  at 

thi 

the  output  of  the  i  element  will  be 


p  <k)  =  ei[WmkAT-iA0] 

1  > 

■3 


(7b) 


...  i  [u>  kAT-iA0-IT/2] 
qi  (k)  =  e  J  m 

=  -j  Pi  no 


Denoting  w.  and  w.  „  as  the  in  phase  and  quadrature  coef- 

1,1  1,2 

ficients ,  the  beamformer  output  signal  is 


z  (k)  =  l  w.  p.  (k)  +  w.  q .  (k) 

i  =  0  1 • 1  x  i , z  i 


jtu  kAT  ,  .  .  -jiA0  ,0. 

=  e  m  ).  (w.  i  ~  iw .  )e  (8) 

i  —  n  1  > 1  1  •  z 


Thus,  the  transfer  function  T(AO)  of  the  system  depends  only 
on  the  progressive  phase  shift  A0  of  the  incident  waveform, 


T ( A0)  = 


(Wi_i-jwii2) 


-jiAO 


and  is  independent  of  the  frequency  of  the  received  signal. 

The  processor  shown  in  Fig.  12  is  preferred  over  that 
shown  in  Fig.  15  because  the  additional  degrees  of  free  am 
--  e.g.,  the  filter  coefficients  --  can  be  utilized  to  ob¬ 
tain  greater  noise  and  interference  rejection.  Although  fre¬ 
quency  selectivity  within  the  receiver  bandwidth  for  signals 
received  from  the  desired  direction  is  clearly  not  a  desirable 
characteristic,  such  selectivity  in  other  directions  allows 
deeper  nulls  in  both  spatial  and  temporal  domains  and  con¬ 
sequently,  decreased  output  noise  power.  Examples  of  this 
phenomenon  are  presented  in  the  sections  following.  It 
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and 


[5] 

should  be  noted  that  constrained  digital  processors 

r  g  I 

associated  adaptive  algorithms  have  been  presented  which 
force  the  system  shown  in  Fig.  12  to  provide  a  flat  frequency 
response  in  the  desired  direction.  These  methods,  however, 
require  a  greater  number  of  digital  calculations  per  adaptive 
step  than  the  unconstrained  procedure  used  in  both  the  pre¬ 
sent  and  the  previous  study^^.  The  constraint  imposed  on 
the  system  also  causes  a  corresponding  reduction  in  null 
depth  in  other  directions.  A  complete  comparison  of  constrained 
and  unconstrained  methods  using  experimental  data,  however,  was 
not  included  within  the  scope  of  the  present  program. 


B  NON  ADAPTIVE  DETERMINATION  OF  BEAMFORMER  COEFFICIENTS 

It  is  presumed  in  this  research  that  the  subarrays  shown 

in  Fig.  12  have  been  steered  in  the  direction  of  the  desired 

target  return  to  the  extent  that  this  signal  appears  in  phase 

at  the  inputs  to  the  digital  filters.  Thus,  the  kth  sample 
th 

from  the  i  subarray  may  be  expressed  as  the  sum  of  a  target 
signal  term  S (k)  and  an  unwanted  component  n. (k) 

xi  (k)  =  S(k)  +  n± (k)  .  (10) 

Note  that  n^(k)  is  the  result  of  noise,  interference,  and 
clutter  returns  received  by  the  ith  subarray.  The  objective 
of  the  beamforming  system  is  to  produce  an  output  signal 
which  closely  approximates  S(k)  —  i.e.,  an  output  which  does 
not  contain  components  produced  by  the  n±{k)  .  The  conven¬ 
tional  approach  to  beamforming  utilizes  only  the  first  set  of 

weights  shown  in  Fig.  12.  Thus,  ail  other  weights  have  zero 
value , 

Wi,£  =  0  1  =  0,1,  ..  K-l;  £  =  1,2,  ..L-l 

(ID 

and  the  first  set  form  a  Tchebychev  or  other  desired  amplitude 
taper.  For  example,  in  the  four  subarray  system  of  interest 
in  the  present  study ,  the  Tchebychev  weights  are  given  by 


w2,0  =  1.0  (12) 
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An  alternative  method  of  beamforming  is  the  array 
processor  studied  in  this  research  in  which  the  filter  co¬ 
efficients  are  chosen  to  minimize  the  average  squared  dif¬ 
ference  between  the  beamformer  output  signal  z(k)  and  a 
desired  output  signal  d(k).  The  desired  response  is  taken 
as  the  target  signal  at  a  delay  of  one-half  of  the  digital 
filter  length  used  to  process  each  subarray  output.  Thus, 

d(k)  =  S (k-L/2) .  (13) 

_2 

and  the  quantity  to  be  minimized  is  given  by  e  ,  where 

e  2  =  E {  (d(k)-z(k)  ]2}.  (14) 


The  set  of  filter  coefficients  which  minimizes  (14)  is 
derived  in  Ref.  1  and  may  be  expressed  in  matrix  notation  as 


W*=  Pxd 


(15) 


where  W*is  the  vector  of  optimum  filter  coefficients 


* 

w  = 


w 


0,0 


w 


1/0 


w 


K  •  0 


w 


K,L 


(16) 


The  square  matrix  P  is  the  matrix  of  cross-correlation 

XX 


values  of  the  data  in  the  processor  and  P  ^  is  the  steering 
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vector  formed  as  the  average  of  the  product  of  the  desired 
array  response  and  the  data  vector  in  the  processor.  By 
assuming  the  S(k)  and  n^(k)  are  uncorrelated  in  (10),  the 
vector  Pxd  is  completely  determined  from  a  knowledge  of  the 
spectral  properties  of  the  target  signal.  For  example,  if 
S(k)  is  sinusoidal  at  frequency  radians/sec  and  the  ADC 

sampling  interval  is  AT  seconds, 


cos  (ui^AT  L/2) 
cos  (uodAT  L/2) 


cos  (u)dAT  L/2) 
cos  [u3dAT(L-l)  /2] 


(17) 


cos  (-w  AT  L/2) 


In  summary,  the  optimum  processor  (16)  may  be  calculated 
from  a  knowledge  of  the  received  data  cross  correlation  values 
and  the  spectral  properties  of  the  desired  beamformer  output. 

This  method  of  selecting  the  coefficient  values  in  Fig.  12  will 
result  in  a  combined  spatial/temporal  filter.  The  precise 
nature  of  this  filter  is  determined  by  the  properties  of  the 
n^(k)  components  in  (10).  In  general,  however,  the  processor 
will  have  a  broadside  mainlobe  at  frequency  uid  while  placing 
pattern  nulls  in  the  direction  and  frequency  of  the  stronger 
interfering  noise  sources.  This  is  in  contrast  with  the 
conventional  beamformer  defined  in  (11)  and  (12)  which  has 
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broadband  nulls  formed  in  predetermined  directions,  inde¬ 
pendent  of  the  spatial  nature  of  the  interference. 


C  ADAPTIVE  DETERMINATION  OF  BEAMFORMER  COEFFICIENTS 


For  a  given  desired  frequency  of  the  target  signal,  the 

optimum  set  of  filter  coefficients  given  in  (15)  is  completely 

determined  from  a  knowledge  of  the  cross-correlation  matrix 

R  .  The  elements  of  this  matrix  are  unknown  a  priori  but 
xx  ~  - 

may  be  estimated  directly  using  finite  time  averages  on  the 


received  signals.  For  example,  an  estimate  of  the  cross¬ 
correlation  between  the  kth  sample  from  the  ith  subarray 
x. (k)  and  the  (k-2)nd  sample  from  the  j  subarray  may  be 


computed  as 


(18) 


There  are  several  difficulties  associated  with  this  approach. 
First,  the  estimate  (18)  is  the  average  over  the  previous  N 
samples  of  received  data,  a  value  chosen  sufficiently  large 
to  provide  a  good  statistical  estimate  but  not  so  large  as 
to  mask  the  effects  of  slowly  varying  noise  parameters  such 
as  direction  of  arrival  of  the  interference.  Selecting  the 
appropriate  value  is  not  an  easy  task.  A  second  difficulty 
arises  in  the  use  of  (18)  to  track  time  varying  intereference . 
Unless  a  large  amount  of  data  memory  is  available,  the 
estimates  (18)  must  be  applied  to  data  received  after  the 
estimates  are  formed.  Thus,  after  each  block  of  N  samples 
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is  used  to  compute  cross-correlation  values,  the  optimum 
coefficients  (15)  are  calculated  and  applied  to  process 
succeeding  data  samples.  At  some  later  time,  new  correla¬ 
tion  values  are  computed  and  the  filter  coefficients  are 
reset  to  the  appropriate  values.  Since  the  filter  is  always 
applied  to  data  at  a  time  after  that  used  to  compute  cor¬ 
relation  values,  it  is  necessarily  suboptimum  for  the  case 
of  time-varying  data. 

The  adaptive  method  of  computing  filter  coefficients 
used  in  this  and  the  previous  study  overcomes  these  dif¬ 
ficulties  in  that  the  coefficients  are  updated  as  each  new 
data  sample  is  received.  A  single,  easily  adjusted  parameter 
in  the  algorithm  controls  the  rate  of  adaptation  or,  equiva¬ 
lently,  the  correlation  averaging  time  of  the  processor.  The 
algorithm^  is 


W  (k+1)  =  W(k)  +  U[Pxd-  z  (k)  X  (k)  ]  (19) 

where  W(k)  is  the  set  of  filter  coefficients  used  to  produce 
th 

the  k  sample  of  the  beamformer  output  signal  z(k)  and  X(k) 
is  the  set  of  data  samples  observed  at  these  coefficients, 
i.e.  ,  f—  -j 


xo(k) 

x1(k) 


X  (k)  = 


x„  ( k ) 

xQ (k-1) 


(20) 


xR (k-L+1) 


The  derivation  of  this  algorithm  is  based  on  a  modi¬ 
fication  of  gradient  descent  techniques.  Details  of  the 
derivation  may  be  found  in  references  1  and  7.  Repeated  use 
of  (19)  results  in  a  sequence  of  coefficient  values  W(k) , 
W(k+1) ,  ...  which  may  be  characterized  as  random  vectors. 

For  the  case  of  stationary  input  statistics,  the  mean  value 
M^(k+1)  of  the  coefficient  vector  at  time  k+1  may  be  computed 
from  (19)  as 


M^k+1)  =  E  [W  (k+1)  ] 

=  E[W(k)  ]+y{Pxd-E[z(k)X(k)  ]} 

-  Vk>+“tPxd-RxxMW(k)>-  (21) 

Evaluation  of  the  last  term  in  (21), 

E  [z  (k)  X  (k)  ]  =  E  [x(k)  XT  (k)  W  (k)  ] 

=  E  [X  (k)  XT  (k)  ]  E  [W  (k)  ] 

=  RxxVk)'  (22) 

is  based  on  the  assumption  that  X(k)  is  a  stationary,  zero- 
mean  gaussian  vector  and  that  the  resulting  coefficients 
are  also  gaussian  [8] .  Further  detailed  discussion  of  these 
assumption  may  be  found  in  references  [9]  and  [10] . 

The  difference  equation  (21)  for  the  mean  coefficient 
vector  may  be  expressed  in  the  form 

Vk+1)  =  tI-^Rxx1Mw(k)+1JPx.r  <23) 
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Repetitive  application  of  (23)  using  an  initial  mean  value 
N^tO)  yields 

M^k)  =  W*+[I-MRxx]kMw(0)  .  (24 

Thus,  the  mean  value  of  the  adapted  coefficients  obtained 
by  iteration  converges  to  the  optimum  set,  provided  that  the 
term  in  square  brackets  converges  to  zero.  It  can  be  shown 
[9]  that  such  convergence  is  guaranteed  provided  that 

0  <  II  <  t-2-  -  (2’ 

max 

where  A  is  the  maximum  eigenvalue  of  the  autocorrelation 
max 

matrix  R  •  A  more  practical  bound  may  be  obtained  by 

XX 

utilizing  the  positive  definite  properties  of  this  matrix. 


Thus 


0  <  U  <- 


K-L-C 


where  a2  is  the  maximum  variance  observed  in  the  array  signals, 
X 

also  insures  convergence. 

For  the  purposes  of  this  study,  a  digital  automatic  gain 
control  was  employed  at  each  element  output  to  ensure  unit 
variance  for  all  channels.  The  proportionality  constant  was 


selected  as 


K-L  f 


where  KL  is  the  total  number  of  adaptive  coefficients  in  the 
beamformer  and  0<a<2  to  ensure  convergence. 

A  measure  of  the  adaptive  time  constant  for  this  processor 
can  be  obtained  from  (24) .  In  general,  the  adaptation  transient 
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will  consist  of  a  sum  of  terms  of  the  form 


(l-P^)*  (28) 

th 

where  0<A.<A  is  the  i  eigenvalue  of  R  .  The  adaptive 

1  luclX  XX 

time  constant  .  is  defined  as  that  value  of  k  which  reduces 

cL 

the  transient  term  to  1/e  of  its  original  value.  Thus,  from 
( 28) ,  we  obtain 


_ -1  <  T  <  .  -1 _ 

ln(l-yA  )  —  a  —  ln(l-yA  .  )’ 

max  min 


(29) 


These  bounds  may  be  simplified  using  (27)  and  the  fact  that 

the  minimum  eigenvalue  of  the  autocorrelation  matrix  cannot 

2 

be  less  than  the  total  power  aN  of  the  uncorrelated  noise  com¬ 
ponent  at  the  processor  inputs.  Further,  the  maximum  eigen¬ 
value  cannot  exceed  the  imaximum  variance  times  the  number  of 


filter  coefficients  and  when  a  -  1,  the  bounds  (29)  are  given 

X 


by 


-1 


-1 


ln(l-a) 


(30) 


It  has  been  observed  in  practice  that  a  reasonable  single 
estimate  for  A^  is  obtained  by  using  an  "average"  eigenvalue. 
Since  all  eigenvalues  must  sum  to  the  trace  of  R  [7] 

XX 


l  A.  =  KL(T  =  KL,  (31) 

i=l  1  x 


and  an  average  value  is  A  =1.0  when  o 
^  av 

responding  time  constant  is  then  given  by 


2 

x 


1.0.  The  cor- 


T 

a 


-1 


inU-  ^  ) 


(32) 
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For  example,  a  processor  with  five  adaptive  weights  at  each 
of  the  four  subarrav  outputs  can  be  expected  to  yield  a  time 
constant  of  about  200  samples  when  a  =  0.10  or  1000  samples 
when  a  =o.02.  At  a  2.0  kHz  sampling  rate,  these  times  cor¬ 
respond  to  0.1  and  0.5  seconds,  respectively.  It  should  be 
noted,  however,  that  rates  considerably  faster  have  beei 

observed  in  practice,  particularly  for  the  case  of  highly 
coherent  —  e.g.,  sinusoidal  --  signals.  Examples  of  this 
phenomenon  are  presented  in  the  sections  following. 

The  discussion  presented  above  has  shown  that  the 
mean  values  of  the  coefficients  obtained  using  the  adaptive 
algorithm  (19)  will  converge  to  the  optimum  values  W*(10) 
provided  that  two  conditions  are  met 

i)  the  proportionality  constant  satisfies  (26) 
ii)  the  vector  used  in  (19)  is  identical  to  the  true 
cross  correlation  vector  defined  in  (17). 

It  can  be  further  shown  [9] -[10],  that  the  variance  of  the 
adapted  solution  can  be  made  arbitrarily  small  through 
sufficiently  small  choice  of  y •  Of  course,  reduced  values 
of  y  imply  corresponding  increases  in  the  convergence  rate 
as  shown  oy  Eq.  (29) . 

D  ADAPTATION  FOR  FM/CW  WAVEFORMS 

The  adaptive  algorithm  described  above  was  designed  for 
operation  in  a  CW  environment  and  as  such  can  be  expected  to 
perform  well  for  both  the  forward  scatter  and  backscatter  CW 
data  files  recorded  during  this  study.  The  data  format  of 
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interest  in  an  operational  HF  backscatter  radar,  however, 
is  that  of  FM/CW  transmissions.  Specifically,  the  trans¬ 
mitted  signal  consists  of  a  sinusoid  in  which  the  instan¬ 
taneous  frequency  is  changed  linearly  with  time  for  the  dura¬ 
tion  of  one  sweep.  At  the  completion  of  each  sweep,  the  fre¬ 
quency  is  returned  to  its  initial  value  with  a  phase  coin¬ 
cident  with  that  used  at  the  beginning  of  the  previous  sweep. 
Figure  5  illustrates  the  instantaneous  time  history  of  the 
specific  FM/CW  format  used  in  this  study.  The  sweep  duration 
of  1/60  sec  provided  32  digital  samples  per  sweep. 

The  receiver  output  signal  generated  by  a  stationary 
target  at  range  rd  consists  of  a  constant- frequency  sinusoid 
which  may  be  set  to  a  prescribed  value  fd  through  appropriate 
choice  of  the  local  oscillator  waveform.  An  example  is 
presented  in  Fig.  16.  Note  that  for  a  brief  portion  of  the 
sweep  ( t d  sec.)  the  beat  frequency  falls  outside  of  the 
receiver  passband .  The  time  waveform  observed  at  the  receiver 
output  for  this  signal  is  shown  in  Fig.  17.  Dashed  lines 
are  used  to  indicate  portions  of  the  sweep  in  which  the 
frequency  is  outside  the  passband  of  the  receiver.  For  the 
parameters  used  in  this  study,  1kHz  receiver  bandwidth  and 
8  1/3  kHz  sweep  range,  the  maximum  allowable  time  period 
outside  the  passband  is  1/(8  1/3)  of  the  sweep  period  or 
about  2  milliseconds. 

Inspection  of  17  shows  that  the  receiver  output  waveform 
for  a  stationary  target  consists  of  a  sinusoid  during  each 
sweep  with  the  phase  of  the  sinusoid  held  constant  with 
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Fig.  17  Receiver  output  time  waveform  for 
stationary  target  shown  in  Fig.  16 


respect  to  the  start  time  of  the  sweep.  A  moving  target  at 
the  same  range  would  produce  a  similar  waveform  except  that 
the  phase  reference  of  the  sinusoid  would  drift  with 
respect  to  the  start  time  of  the  sweep.  In  Fig.  18,  succes¬ 
sive  sweeps  are  arranged  vertically  to  illustrate  this  property. 
Note  that  the  above  discussion  is  based  on  the  assumption  that 
the  doppler  frequency  caused  by  target  motion  produces  only 
a  small  frequency  change  during  one  sweep.  This  assumption 
is  consistent  with  the  targets  of  interest  in  the  HF  environ¬ 
ment  . 

The  signal  observed  at  the  output  of  an  FM/CW  radar 
during  any  one  sweep  may  thus  be  approximated  by  a  sinusoid 
at  frequency  f^,  determined  by  the  target  range,  and  at 
phase  <^,  determined  by  the  target  range  and  doppler,  and 
by  the  sweep  number.  The  transient  observed  after  frequency 
flyback  consists  of  both  an  amplitude  change  and  a  phase  shift. 
An  adaptive  system  which  is  designed  to  operate  in  the  pre¬ 
sence  of  a  sinusoidal  desired  output  signal  can  be  expected 
to  perform  well  during  the  mid-sweep  portions  of  the  re¬ 
ceived  signals.  The  effect  of  the  flyback  transient  on  adap¬ 
tive  performance,  however,  is  not  so  easily  analyzed.  One 
method  of  avoiding  the  influence  of  these  transients,  is  to 
cease  adaptation  for  the  period  in  which  the  transients 
reside  in  the  processor.  This  procedure  has  beem  implemented 
and  tested  on  recorded  FM/CW  data.  The  results  indicate  that 
adaptive  performance  is  unaffected  by  these  transients.  Spe¬ 
cific  examples  are  presented  in  the  sections  following. 


Sweep 

number 

1 


2 


3 


a)  stationary  target 


i 


b)  moving  target 


Fig.  13  Receiver  output  waveforms  for  fixed  and 
moving  targets  in  an  FM/CW  system 
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V. 


RESULTS 


The  results  presented  below  are  divided  into  three  cat¬ 
egories  according  to  the  experimental  data  configuration  used 

to  record  the  data  —  i.e.,  forward  scatter  CW  tests,  back- 
scatter  CW  tests  and  backscatter  FM/CW  tests.  In  each  case 
examples  are  presented  comparing  the  performance  of  the 
adaptive  beamformer  with  that  of  the  conventional  beamformer. 
The  results  presented  were  derived  using  a  wide  variety  of 
adaptive  parameters  and  input  signal  configurations.  All 
computations  were  carried  out  using  the  CDC  6400  computing 
facility  at  the  University  of  Colorado. 

A.  FORWARD- SCATTER  CW  TESTS 

In  order  to  demonstrate  the  spatial-temporal  filtering 
capabilities  of  the  tapped-delay-line  processor  shown  in 
Fig.  12,  a  number  of  non-adaptive  computations  were  carried 
out  using  the  optimal  coefficient  relationship  defined  in 
Eq .  (15) .  In  each  case,  a  total  of  128  data  samples  per 

element  were  used  to  estimate  the  autocorrelation  matrix 
values.  All  time  averages  were  carried  out  using  the 
method  shown  in  Eq.  (18) .  The  inverse  of  the  resulting 
matrix  was  then  multiplied  by  the  steering  vector  in  Eq.(17). 
A  desired  output  frequency  wd/2u=  500  Hz  corresponding  to 
the  downcon ve rted  CW  Bearden  signal  was  used  in  these  cal¬ 
culations  . 

The  input  data  were  taken  from  File  3  and  consisted  of 
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a  low  SNR  forward  scatter  CW  signal  received  in  the  pre¬ 
sence  of  teletype  interference  (see  Table  I,  Section  III) . 
Figure  8  shows  the  azimuth/f requency  content  of  the  first 
128  sample  data  block  recorded  on  this  file.  Optimal  filter 
coefficient  values  were  computed  for  processors  containing 
1,  3,  5,  and  9  taps  per  array  element.  A  display  of  the 
azimuthal  response  of  the  resulting  processor,  plotted  as  a 
function  of  frequency  within  the  received  bandwidth,  was 
generated  using  two-dimensional  Fourier  Transform  techniques. 

The  calculation  proceeds  by  arranging  the  filter  co¬ 
efficients  in  a  two-dimensional  array  with  each  row  con¬ 
taining  the  coefficient  values  for  the  tapped-delay-line 
attached  to  a  particular  antenna  element.  Thus,  with  ref¬ 
erence  to  Figure  12  and  Eq .  (16) ,  the  matrix  is 


Wq  o  WQ  1  ••  Wq 

Wl,0  W1,L-1 


(33) 


WK-1 , 0  WK-1,L-1  , 

Transforming  the  rows  and  columns  of  (33)  leads  to  a  matrix 
in  which  each  row  is  the  frequency  response  of  the  processor 
in  a  particular  azimuthal  direction.  The  azimuthal  spacing 
is  determined  by  the  original  received  frequency  and  the 
subarray  separation.  It  should  be  noted  that  the  validity 
of  this  procedure  is  based  on  the  assumption  that  the  sub- 
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array  spacing  in  wavelengths  remains  constant  over  the 
bandwidth  of  the  received  signals.  Thus,  the  technique  may 
readily  be  applied  to  narrowband  HF  arrays  operating  at 

bandwidths  less  than  50  kHz. 

The  azimuth/frequency  responses  obtained  with  this 
method  for  the  optimal  filters  derived  from  the  first  128 
data  samples  in  File  3  is  shown  m  Figs.  19-22.  In  each  case, 
the  product  of  the  data  transform  in  Fig.  8  and  the  appro¬ 
priate  processor  transform  is  also  illustrated.  The  desired 
CW  signal  is  indicated  by  the  vertical  arrow.  Inspection 
of  this  series  of  plots  shows  that  the  received  signal-to- 


noise  ratio,  determined  by  a  summation  over  the  product 


pattern,  improves  as  the  number  of  filter  coefficients  in¬ 
creases.  The  improvement  results  from  an  increased  selecti¬ 
vity  in  both  azimuthal  and  temporal  filtering.  Thus,  a 
comparison  of  Figure  19  with  Figure  22  shows  that  the  strong 
teletype  signal  is  eliminated  at  all  azimuths  in  the  latter 
case  but  only  on  two  azimuths  in  Fig.  19. 


An  indication  of  the  time  variability  of  the  optimal 


pattern  for  the  5  tap  per  element  case  is  shown  in  Fig.  23. 

In  this  series,  optimal  coefficients  were  computed  for  4 
successive  blocks  of  128  data  samples  from  File  3.  Sub¬ 
stantial  plot-to-plot  variation  is  evident,  particularly 
in  the  sidelobe  structure.  For  example,  compare  the  frequency 
response  in  the  azimuthal  direction  for  Figs.  (23b)  and  (23c). 
This  variation  cannot  be  attributed  to  statistical  fluctua- 
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tions  in  the  autocorrelation  estimates  and  is  therefore 
indicative  of  changes  in  the  azimuth/ temporal  structure  of 
the  received  data. 

The  P-vector  adaptation  algorithm  (19)  was  also  imple¬ 
mented  on  these  data.  At  intervals  of  128  samples,  the 
adapted  weight  vector  was  read  out  and  azimuth/frequency 
displays  were  generated.  Figure  24  illustrates  the  results 
obtained  at  four  successive  intervals  by  adapting  on  the 
same  data  as  used  to  compute  the  optimal  patterns  in  Fig.  23. 
The  adapted  processor  contained  5  coefficients  per  element 
and  the  value  of  a  selected  [see  Eq.  (27)]  was  0.4.  Al¬ 
though  the  adapted  and  time-averaged  plots  are  markedly 
similar,  differences  in  the  sidelobe  structure  can  be  de¬ 
tected.  These  differences  can  be  attributed  to  the  fact 
that  time-averaging  is  a  block-time  process  taken  over  128 
samples  and  adaptation  is  equivalent  to  exponential  averaging 
with  the  length  of  the  average  being  determined  by  the 
adaptive  time  constant.  Using  the  method  presented  in  Section 
IV,  the  average  time  constant  for  this  set  of  parameters 
is  about  50  samples. 

Time  history  plots  of  the  received  signal-to-noise  ratio 
for  File  3  were  computed  for  both  conventional  and  adaptive 
processors.  The  SNR  was  computed  from  time  averages  on  the 
received  data  as  the  ratio  of  output  sinusoidal  power  at  the 
desired  frequency  to  total  output  power.  Thus,  the  maximum 
possible  output  SNR  was  unity  or  zero  db.  Figures  25  and  2b 
illustrate  this  quantity  during  the  first  140  m.sec.  of 
File  3  for  adaptive  processors  containing  one  weight  and 
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y  of  output  SNR  for  forward-scatter  CW  test 
ficient  per  element 


nine  weights  per  element,  respectively.  The  large  transient 
observed  near  45  m.  sec  was  produced  by  setting  the  adapted 
weights  to  zero.  In  both  cases,  the  observed  recovery  time 
was  about  20  m.  sec.  A  summary  of  the  output  SNR  improvement 
provided  by  adaptive  processing  on  File  3  is  shown  in  Table 
III.  Gain  only  adaptation  provides  on  the  order  of  3  db  im¬ 
provement  while  nine  coefficients  per  element  produces 
greater  than  10  db .  improvement. 

B  BACKSCATTER  CW  TESTS 

Output  SNR  calculations  were  also  carried  out  for  data 
recorded  using  backscatter  CW  transmissions.  Figures  9  and 
10  illustrate  the  azimuth/frequency  content  of  these  data 
under  both  high  (File  7)  and  low  (File  8)  input  SNR  condi¬ 
tions.  a  desired  signal  was  generated  using  a  100  Hz  offset 
repeater  at  Bearden.  The  equivalent  downconverted  value  at 
the  Los  Banos  receiver  output  was  then  400  Hz.  Note  that 
an  image  frequency  line  appears  at  600  Hz  due  to  the  use 
of  a  double  sideband  modulator.  For  purposes  of  the  follow¬ 
ing  discussion,  only  the  400  Hz  line  was  taken  as  the  desired 
output  signal. 

A  time  history  of  the  SNR  measured  at  the  output  of 
both  a  conventional  beamformer  and  an  adaptive  signal  pro¬ 
cessor  operating  on  File  8  with  5  adaptive  weights  per  an¬ 
tenna  element  u  =  0.4  is  shown  in  Fig.  27.  The  SNR  in 

this  plot  was  computed  as  the  ratio  of  output  power  at  the 
desired  400  Hz  frequency  to  total  output  power.  In  effect, 
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Ill:  Summary  of  output  SNR 
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the  quantity  measured  is  the  ratio  of  signal  to  interference 
plus  clutter  plus  noise  power.  The  adaptive  coefficients 
were  set  to  zero  at  time  t  =  40  m.sec  for  purposes  of  mea¬ 
suring  the  adaptive  recovery  time.  A  summary  of  the  improve¬ 
ment  offered  by  the  adaptive  processor  with  respect  to  con¬ 
ventional  beamforming  is  presented  in  Table  IV. 

The  results  obtained  in  the  presence  of  CW  backscatter 
transmissions  were  markedly  similar  to  those  measured  for 
the  forward  scatter  case.  Larger  improvements  in  the  back¬ 
scatter  case  are  primarily  due  to  the  clutter  rejection  pro¬ 
vided  by  the  adaptive  processor.  Figure  28  illustrates 
this  rejection  in  more  detail.  A  sequence  of  four  azimuth/ 
frequence  response  curves  were  computed  using  the  adapted 
filter  coefficients  obtained  during  the  experiment  shown 
in  Fig.  27.  The  response  at  the  clutter  frequency  (500  Hz) 
is  suppressed  at  all  azimuths  except  borsite. 

These  experiments  were  carried  out  to  investigate  the 
performance  of  the  adaptive  processor  in  the  presence  of  an 
undesired  signal  (ciutter)  which  was  highly  correlated  with 
the  desired  signal.  Continuous  wave  transmissions  were  used 
for  purposes  of  providing  simple  signal  formats  and  were  not 
intended  to  be  representative  of  operational  systems.  The 
results  indicated  that  adaptive  processing  works  well  in 
the  presence  of  clutter  and  no  new  operational  difficulties 
were  encountered. 
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EV :  Summary  of  output 


Sequential  adaptive  processor  response  plots  derived  from 
File  8  with  five  coefficients  per  element. 


C  FM/CW  RADAR  EXPERIMENTS 


The  data  recorded  in  Files  10,  11,  and  12  were  taken 
using  the  FM/CW  transmission  format  shown  in  Fig.  5.  These 
data  were  processed  using  both  a  conventional  beamformer  and 
the  adaptive  configuration  shown  in  Fig.  12.  Range/doppler 
maps  were  computed  from  the  processed  outputs  using  two- 
dimensional  Fourier  Transform  techniques.  Specifically, 
denoting  the  kth  sample  of  a  processed  beam  output  as  z(k) , 
a  two-dimensional  matrix  is  constructed  in  which  the  nth  row 
contains  the  beam  output  samples  during  the  nth  sweep.  With 
a  60  Hz  sweep  rate  and  a  sampling  frequency  of  1.92  kHz, 
there  are  precisely  32  samples  per  sweep  and  the  matrix  is 
given  by, 

z (0)  z ( 1 )  .  .  z ( 31) 

z  ( 32)  .  .  .  z  (63) 


(34) 


z [ (N- 1) 32] 


z  ( 32N-1) 


The  number  of  sweeps  N  determines  the  dwell  time  for  the  map. 
In  the  examples  shown  below,  a  value  N  =  32  is  used  providing 
a  dwell  time  of  0.533  seconds. 

Transforming  the  rows  and  columns  of  this  matrix  pro¬ 
duces  a  matrix  of  complex  numbers  with  dimension  N  x  16 . 

(The  reduction  from  32  to  16  is  due  to  the  use  of  real-valued 
samples  in  the  initial  array.)  After  transforming,  the  mag- 
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nitude  of  the  values  in  the  mth  column  of  the  resulting 
matrix  may  be  interpreted  as  the  doppler  response  in  the  mth 
range  bin.  Thus,  with  N  =  32,  the  process  produces  a  range/ 
doppler  map  with  16  range  bins  and  32  unambiguous  doppler 
cells.  For  the  parameters  used  in  this  experiment,  the  range 
bin  separation  was  18.75  km  and  the  doppler  cells  were  sep¬ 
arated  by  1.875  Hz. 

A  range/doppler  display  was  generated  from  the  magnitudes 
of  the  transformed  matrix  by  normalizing  with  respect  to 
the  root-mean-square  noise  floor.  Figure  29  shows  the  re¬ 
sults  of  five  successive  dwells  computed  using  the  data  in 
File  11.  In  each  case,  maps  derived  from  both  conventional 
and  from  adaptive  beamforming  outputs  were  computed.  The 
two  outputs  were  generated  simultaneously  using  identical 
input  signals.  The  adaptive  processor  consisted  of  five 
coefficients  per  element  and  a  value  of  a  =  0.1  The  steer¬ 
ing  vector  frequency  used  was  500  Hz  and  corresponds  to  the 
range  indicated  by  the  arrow.  These  data  contain  a  clear, 
interference-free  repeater  signal  located  in  the  8th  range 
cell.  The  target  doppler  appears  at  both  +16.9  and  -16.9 
Hz  due  to  the  use  of  double  sideband  modulation  in  the  re¬ 
peater.  Ground  clutter  appears  in  three  doppler  cells  and 
is  clipped  at  a  value  ten  times  greater  than  the  RMS  noise 
floor.  Display  amplitudes  are  presented  on  a  linear  scale 
and  Hamming  windows  were  applied  to  the  data  in  both  dimen¬ 
sions  prior  to  transforming. 

A  comparison  of  adaptive  and  conventional  processing 
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+  30  Hz. 


15 


Fig.  29  a)  Range/doppler  displays  derived  from  File  11, 
Dwell  #  1. 
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b)  Range/doppler  displays  derived  from  File  11 
Dwell  #  2. 


DOPPLER 


ADAPTIVE 


Fig.  29  d)  Range/doppler  displays  derived  from  File  11, 
Dwell  #  4. 


ADAPTIVE 


29  e)  Range/doppler  displays  derived  from  File  11, 
Dwell  #  5. 


using  data  with  a  lower  signal-to-noise  ratio  is  shown  in 
Fig.  30.  Data  for  these  results  were  obtained  from  File  10. 

The  filter  parameters  used  to  generate  this  display  were 
identical  to  those  used  in  Fig.  29.  The  five  dwells  shown 
are  time  sequential  with  adaptive  and  conventional  outputs 
computed  simultaneously. 

As  discussed  in  Section  IV  D,  the  phase  and  amplitude 
transient  which  occurs  during  frequency  flyback  in  the  FM/ 

CW  transmission  may  be  expected  to  affect  the  performance 
of  the  adaptive  processor.  Figure  31  presents  results  ob¬ 
tained  with  a  processor,  in  which  adaptation  occurred  only 
during  the  middle  portions  of  the  frequency  ramp  --  i.e. 
during  those  periods  in  which  the  transient  was  absent.  The 
two  range/doppler  maps  shown  were  computed  using  the  same 
data  and  filter  parameters  as  used  for  the  first  dwells  shown 
in  Fig.  30.  No  apparent  advantage  was  detected  for  the  method. 

The  effect  of  variation  in  the  adaptive  time  constant 
is  illustrated  in  Figs.  32-39.  This  sequence  of  eight 
successive  dwells  was  computed  using  the  pseudo-operational 
data  recorded  in  File  12.  Examples  of  adaptation  using  a 
5  weight  per  element  processor  with  a  =  0.02,  0.10  and 
0.50  are  presented.  The  steering  vector  frequency  selected 
was  340  Hz  aid,  in  effect,  optimizes  the  processor  for  all 
doppler  frequencies  at  the  range  point  indicated  by  the  arrow. 
These  results  indicate  that  wide  variations  in  the  value  of 
u  produce  markedly  similar  adapted  range/doppler  maps.  The 
average  time  constants  computed  from  (32)  are 
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Fig.  30  b)  Range/doppler  displays  derived  from  File  10, 
Dwell  #  2. 
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CONVENTIONAL 


30  c 


)  Ranae/doppler  displays  derived  from  File  10, 
Dwell  #  3. 


CONVENTIONAL 


ADAPTIVE 


d)  Range/doppler  displays  derived  from  File  10, 
Dwell  #  4. 


30  e)  Range/doppler  displays  derived  from  File  10, 
Dwell  #  5. 


rom 


Range/doppler  displays  derived  from  File  12,  Dwell 


RANGE 


Adaptive  processor. 


t  value 


t  (m. sec) 

0.02  520 

0.10  103 

0.50  20 

'  1US '  the  smallest  value  of  a.  used,  the  time  constant 
was  approximately  equal  to  the  dwell  period.  Rapid  adaptation 
nay  induce  an  output  noise  component  which  is  produced  by 
fluctuation  of  the  processor  coefficients.  An  example  of 
this  effect  at  a  =  0.5  is  evident  in  Fi  (34c). 

Cor.pu tations  were  also  carried  out  to  determine  the 
decree  to  which  the  coefficient  values  were  oscillating  due 
to  adaptation  durincj  this  experiment.  Mean  and  root-mean- 

square  (RMS)  values  for  the  coefficients  were  derived  using 
a  128  sample  average.  The  results  obtained  with  a  =0.1 
are  sunmarized  in  Table  V.  Similar  tables  were  constructed 
for  u  =0.5  and  t  =  0.02  and  the  data  for  all  three  a 
vi lues  wore  used  to  compile  an  average  "coefficient  variation" 
v  as  follows 
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where  K  = 

L  = 

M.  . 

0  . 

1 J 

The  values 


no.  of  elements 
no.  of  coefficients/element 
=  mean  of  jth  coefficient  at  ith  element 
=  RMS  of  jth  coefficient  at  ith  element, 
obtained  as  a  function  of  n  and  expressed 


as 


a 


percentage  were 


t  vaLue 
0.0  2 


Y 

2.26? 


88 


0.10 


3.69 


0.50  7.75 

These  results  indicate  that  larger  values  of  a  produce 
a  creator  percentage  coefficient  oscillation,  as  expected 
from  the  known  theory  of  convergence  for  the  P-vector  algorithm. 
Observation  of  the  r ange/dopple r  maps  shown  in  Figs.  32-39 
further  indicates  that  y  values  the  order  of  8%  induce 
detectable  changes  in  the  map.  With  a  =  0.1,  the  percentage 
change  in  coefficients  was  reduced  to  less  than  4%  and  no 
adaptation  effects  were  observed  in  the  ranae/doppler  outputs. 
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VII  DISCUSSION  AND  CONCLUSIONS 


This  report  has  investigated  the  applicability  of  adaptive 
Le conforming  techniques  to  HF  radar  arrays.  It  has  been  shown 
that  a  simple  adaptive  algorithm  can  be  used  to  combine  individual 
received  subarray  signals  within  a  large  array  so  as  to  provide 
an  effective  increase  in  the  received  signal-to-noise  ratio 
relative  to  conventional  beamforming  methods.  The  adaptive 
processor  has  been  shown  to  provide  combined  azimu‘.hal  and  fre¬ 
quency  selectivity  within  a  one  kHz.  bandwidth  centered  at  about 
eleven  MHz. 

Specific  results  were  obtained  for  the  case  of  forward 
scatter  one-hop  ionospherically  propagated  signals  received  in 
the  presence  of  teletype  and  other  interference  and  noise  wave¬ 
forms.  Signal-to-noise  ratio  improvements  observed  for  this 
class  of  signals  varied  between  3dB,  for  the  case  of  one  adaptive 
coefficient  per  subarray  element,  to  about  11  dB  for  the  nine 
coefficients  per  element  case.  Total  time  taken  to  adapt  from 
zero  weight  values  to  the  optimal  solution  was  less  than  40 
m.sec.  for  all  tests. 

Adaptive  experiments  were  also  carried  out  in  the  presence 
of  CW  backscatter  signals.  These  tests  demonstrated  that 
adaptation  can  significantly  increase  the  received  SNR  of  a 
repeater  signal  offset  by  100  Hz  from  the  clutter.  Adaptive 
performance  in  the  presence  of  clutter  was  for  all  practical 
purposes  identical  with  that  observed  in  the  forward-scatter 

experiments.  Greater  than  25  dB  improvements  against  clutter 
and  noise  were  observed. 


Perhaps  the  most  significant  results  were  obtained  using 
FM/CW  sweep  frequency  transmitter  waveforms.  Range/doppler  maps 
were  computed  from  receiver  array  output  signals  using  both 
conventional  and  adaptive  beamforming  methods  to  combine  the 
subarray  outputs.  Tests  were  conducted  using  both  known  repeater 
signals  and  using  targets  of  opportunity.  In  all  cases,  the 
range/doppler  display  produced  with  adaptive  beamforming  was 
significantly  improved  over  that  generated  by  conventional 
steering.  Specific  improvements  in  excess  of  12  dB  within  a 
selected  range/doppler  cell  were  observed  on  several  occasions 
(Figs .  33  and  34 )  . 

The  adaptive  algorithm,  employed  in  these  studies  was 
originally  designed  for  use  in  passive  sonar  receiving  arrays 
and  has  not  been  optimized  for  use  in  the  FM/CW  radar  environment. 
The  results  obtained,  however,  demonstrate  that  the  procedure  is 
well  suited  to  FM/CW  transmissions  and  steps  taken  in  this 
research  to  further  optimize  the  algorithm  did  not  lead  to 
significant  performance  improvements. 

The  primary  conclusion  drawn  from  the  results  presented 
herein  is  that  time-domain  adaptation  can  significantly  improve 
the  performance  of  an  HF  FM/CW  radar  system.  The  procedure  may 
be  readily  implemented  in  real  time  using  either  special  purpose 
hardware  of  a  fast  general-purpose  computer.  Further  research, 
however,  is  required  to  determine  the  optimal  implementation 
configuration.  This  research  should  consider  both  analog/digital 
hybrid  systems  and  the  all  digital  approach  used  in  the  present 

study.  Additional  investigations  should  also  be  carried  out  to 
determine  the  optimal  subarray  geometry  for  the  system. 
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